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Abstract. We show that spin and fermion representations for solvable quantum 
chains lead in general to different reduced density matrices if the subsystem is not 
singly connected. We study the effect for two sites in XX and XY chains as well 
as for sublattices in XX and transverse Ising chains. 



1. Introduction 

The investigation of entanglement features in many-body quantum systems has been 
a topic of intense research in recent years pQ. In such studies, one divides the system 
in two parts in space and asks how these are coupled in the quantum state. This can 
be answered from the reduced density matrix (RDM) for one of the subsystems. For 
a number of solvable lattice models, these density matrices can be found exactly [2] 
and their spectra determine the entanglement entropy which is a simple measure of 
the coupling. In the ground state of typical systems with short-range interactions, it 
is connected with the interface between the two parts and proportional to its extent, 
although there are corrections to this "area law" for fermions O 0] . 

Commonly, one divides the system into two parts which are singly connected, 
and most results were obtained for this geometry. However, there are a number of 
studies which treat multiply connected subsystems. These were investigated in one 
dimension for free fermions and bosons within a continuum approach [5j [6j and for 
conformally invariant systems [U [7]. On a lattice, spin chains have been treated 
where the subsystem consisted of many single spins separated by a certain number 
of sites [8], of two blocks of spins, again separated OfTO], and of a whole sublattice 
PTj . The latter choice was motivated by the search for phase-transition indicators 
while for two blocks one can use the RDM to study the entanglement between them. 
In [8l [9] the calculation was carried out by transforming the transverse Ising (TI) 
model into fermions and determining the density-matrix eigenvalues from the fcrmionic 
correlation functions [12j [13l HU [15]. The basic finding for isolated spins was that, 
due to the many "interfaces" with the surrounding, the entanglement entropy becomes 
extensive. 

In this note, we want to point out that the reduced density matrices for the 
spin representation and the fermionic representation are, in general, not identical 
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and that correspondingly also the entanglement entropies differ, if one deals with 
disjoint subsystems. This might be surprising at first, but is connected with the non- 
local structure of the Jordan- Wigner transformation. Thus, to obtain a transverse 
spin correlation function, one needs information about a whole string of sites in the 
fermionic picture [16]. This is not necessary if one asks only for fermionic correlations. 
Thus the two RDM's contain different information. In the following section 2, we 
demonstrate this for the simple and analytically solvable case of a subsystem of two 
sites in an XX chain. In section 3 we indicate the generalization to the XY chain. The 
possibly more relevant case of sublattice RDM's is treated in section 4 for XX and TI 
rings and the results are summarized in section 5. 



2. Two sites in an XX chain 



In the following we consider the spin one-half XX chain described by the Hamiltonian 

m 

Here the tr^ are Pauli matrices at site m. In terms of fermionic creation and 
annihilation operators cj„ , c m , H reads [TB] 

H = - J2( c Lc m +i + c] n+1 c m ) (2) 

m 

and corresponds to a simple hopping model. If the system forms a ring, one has to 
take care of a boundary term. In the following we always look at the ground state 
and choose as subsystem the sites 1 and n. 



2.1. Spin RDM 



The RDM for two spins in XX, XY and TI chains has been discussed in many papers, 
see [17l [HI [TH |20j |2~U [221 [23] . It is a 4 x 4 matrix and its elements can be expressed as 
expectation values of proper operators at the two sites, i.e. as two-point correlation 
functions. Since the XX ground state has fixed S% ot = 0, all entries corresponding to a 

change of this value are zero. As a result, the RDM in the basis |++), |H — ), | — +), | ) 

has the simple form 

\ 

(3) 



/ a 



Ps 



V 



a J 



where 2a + 2b = 1 due to tr(ps) = 1- The non-zero matrix elements are given by 

1 



a = p S {++,++) = + + 



Ps(+- +-) = i ((l + af)(l 



-+) = 



where 



a x ± a v . In the ground state, (cr^ 
«=7 + 7W- 



<)), 



and the expressions become 



(4) 
(5) 
(6) 

(7) 
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(8) 

(9) 



In this way, the RDM is expressed completely in terms of standard spin correlation 
functions. This approach has also been used to determine RDM's in the XXZ model, 
see e.g. [53] and references therein. 

The correlation functions can be calculated in the fermionic representation. Then 
one has 

■g( a i a n) = - Ci,nC„,i = _ d n (10) 

where Cij = (c\cj) is the fermionic single-particle function. For the xx correlations 
one finds the determinant [H>] 



(11) 





Cl,2 


Cl : 3 


■ Ci <n 




C2,2 


C*2,3 


■ C*2,n 


2 n-2 


^3,2 


C3,3 


• Cl,n 




C n -l t 2 


Cn-1,3 


Cji — 1 , n 



1/2 and vanish in the ground 



Here the quantities C m ,m are given by C m , 
state, for which the fermionic system is half filled. In two cases the expression reduces 
to a single term, namely if the two sites are nearest neighbours (n — 2), and if they are 
the ends of an open chain. Then the effect of the Jordan- Wigner strings disappears 
and 

1 -{ofo x n ) = C 1 , n (12) 



The "long-distance entanglement" in the latter case has been studied e.g. in [2"Tl I25j . 



2.2. Fermion RDM 

In the fermionic case, one can proceed in exactly the same way. The RDM is 
again a 4 x 4 matrix, this time in the basis specified by the occupation numbers 
|11), |10), 1 1 ) , 1 00} . Since the ground state has fixed total particle number A/", all 
elements corresponding to a change of Af vanish. Therefore p p has again the form §3§ 
and the elements are now given by 

a = (c\ Cl cic n ) = -- Cl n , (13) 

b = {c\c l {l-c\ l c n )) = \+Cl n , (14) 

c = (c{c n ) = Ci,„. (15) 

One sees that a and b are the same as in the spin representation, but c is in general 
different. Instead of the xx spin correlation function, the fermionic one-particle 
correlation function C\ >n appears. This is the basic difference between the two 
representations. Only for the two exceptional cases mentioned above, the expressions 
for c coincide. Then the two spins are nearest neighbours, or can be considered as 
neighbours by bending the open chain to a ring, and correspondingly the subsystem 
is singly connected. 
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The appearance of C\ >n is natural, since with pp one must in turn be able to 
calculate this correlator. That essentially only this quantity enters, could have been 
seen also from the general result for fermionic RDM's [T21 H] 

Pf = -| exp(-W) = - exp(- ^ H ltj c\c 3 ) (16) 

i J — l,n 

where Z ensures tr(pp) = 1 and the matrix Hij in the exponent follows from the 
correlation matrix involving the sites of the subsystem, in our case Ci t \ = C n ^ n = 1/2 
and Ci,„. One can obtain the form ([3]) also from (|16p . but the route taken above is 
much simpler. 



2.3. Entanglement entropies 

/From the structure of p g and p f one reads off the four eigenvalues 

W\ = b + c, W2 = W3 = a, W4 = b — c. (17) 

The entanglement entropy is then given by S = — Wk In Wk ■ Since the correlations 
go to zero for large separations, all Wk approach 1/4 in this limit and S goes to the value 
S = 2 In 2 in both representations. This is the sum of two single-spin contributions 
and also the maximum one can have. 

The complete behaviour of S is obtained by using the correlation function of the 
infinite chain 

a, = sin( t'f 2) (is) 

and calculating the eigenvalues Wk numerically. The result is shown in FigfT] 
According to the previous remarks, the two entropies coincide for n — 2 and for 
n — > oo. In between, the spin entropy Sg always lies below the fermionic one, Sp, 
because the ix-correlations only decay as n -1 / 2 . If the distance n — 1 is even, i.e. 
if the two sites are on the same sublattice, the correlation Ci >n vanishes and Sf 
has the maximum value 2 In 2. Apart from this feature, the asymptotic behaviour is 
determined by c 

S"- 2 In 2 -4c 2 (19) 

which gives a 1/n approach to the limit in the spin picture and a 1/n 2 variation in 
the fermion picture. 



3. Anisotropic chains 



The previous considerations can be generalized to the case of an XY or a TI 
chain. Then the rotational symmetry is absent, but the ground state is built from 
configurations with either an even or an odd number of (+) spins. Then 

(a d\ 

I I ( 20 ) 



V d 



a J 



On reduced density matrices for disjoint subsystems 



5 



<AJ 
00 



1.7 



*. *~ 4 .- + •••+ 
v 



+ fermion 
• spin 




Figure 1. Two-site entanglement in the infinite XX chain as a function of the 
distance n, calculated in the spin and in the fermion basis. 



and the additional matrix element d is given by 

d = Ps(++,—) = l(cria+) (21) 
In the fermionic picture, on the other hand 

d = pp(H,00) = <cJ4) (22) 

Again, both expressions are different, except for n = 2 or for end spins. The presence 
of d splits the two degenerate eigenvalues of p and gives W2,3 = a± d. 

An additional feature arises from the long-range order which exists in the 
anisotropic XY model and in the TI model with strong coupling. In the spin picture, 
the coefficients c and d do not vanish for large separations and therefore the entropy 
does not become the sum of two single-site contributions. If 7 denotes the anisotropy 
of the XY model with zero field, the asymptotic values are [26] 

c = rf=-((7X) = -— (23) 
2 \ 1 n/ 2I + 7 y ' 

and vary between for 7 = and 1/4 for 7=1. Correspondingly, Ss varies between 
2 In 2 and In 2. The fermionic correlation functions, on the other hand, approach zero 
for large separations, as in the XX model, and Sp = 2 In 2 independent of 7. Thus the 
quantities <Sg and Sp have different asymptotic values. Only for end spins in open 
chains, the asympotic values coincide and are in this case related to the surface order. 
This happens also in dimerized XX chains [3TJ [35] . 



4. Sublattice entanglement 

We now turn to the case where the subsystem consists of every second site in a chain. 
For the infinite XX model, it then follows from ([T5]) that the fermionic correlations 
vanish for unequal sites on the same sublattice. This holds also for finite rings, as 
well as for non-homogeneous couplings. Therefore the correlation matrix Cjj on one 
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sublattice is diagonal, Cjj = 1/2 Sij. The RDM obtained via (fH))) is then also 
diagonal, all eigenvalues Wk are equal and Sp = Lin 2 if the system is a ring with 
N = 2L sites. This is the extensivity mentioned in the introduction, and the result has 
a simple interpretation. The ground state |\&) is obtained by filling the single-particle 
eigenstates of ([2]) for momenta \q\ < ir/2. But the corresponding operators can be 
decomposed as 



J - 



1 



ex p H9 2to )4 



1 



L 

N ^ 

m—l 



exp (~iq(2m + 1))<4 



2m+l 



V2 



(4 + 6j) 



Therefore 



I*) 



n 



7i K 



b\)\0) 



(24) 



(25) 



? |<7r/2 

is a product of fermionic triplets formed from modes with the same q on different 
sublattices. Each of them contributes In 2 to the entanglement entropy. 



For the spin RDM, a calculation for N — 4 already shows that the Wk are not all 
equal, but given by Wi = 1/2, W2 = w?, = 1/4, 104 = 0. The corresponding entropy is 
Ss = 3/2 In 2 and thus smaller than Sp — 2 In 2. To investigate the size dependence, 
we have calculated ps numerically by finding the ground state of H in the subspace 



Hot 



for N up to 16 sites. The resulting values are plotted in Fig[2]on the left. 




fermion 
spin 



2.5 - 



N 



14 



Figure 2. Sublattice entanglement entropy for rings of length N = 2L, calculated 
in the spin and in the fermion basis. Left: XX model, right: critical TI model. 
The straight lines represent the respective asymptotic behaviour, see the text. 



The data can be fitted for the larger L by S's(L)/ln2 = s * L + so with 
s = 0.6116(4) and sq = 0.290(5). In addition, there are subleading corrections which 
are different for even and odd values of L/2, related to the slightly alternating be- 
haviour of ps- Thus the extensivity of S holds also in the spin picture, but the value 
is reduced to about 60% of the fermionic one. For open chains, the values for S$ lie 
about 0.5 higher than for the rings, but the slope s is similar. 
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We have also considered the XX chain with enforced dimerization where the 
coupling between sites m and m + 1 is given by J m = 1 + S (— l) m . Then one has 
a finite correlation length if 6 =/= 0. However, the fermionic sublattice entanglement 
is not affected. The expectation values (c\cj) still vanish on the same sublattice and 
one can interpret the result Sf = Lin 2 as before. The spin entanglement, on the 
other hand, is an even function of 8 and depends on the dimerization. In particular 
for 5 = ±1, where the system decomposes into coupled pairs of sites each of which 
contributes In 2, it becomes equal to Sf- 




0.6 1 1 1 

0.5 „ 1 

o 

Figure 3. Sublattice entanglement entropy per site of XX rings as a function of 
the dimerization, calculated in the spin basis for different sizes TV. 

The detailed variation with 5 for fixed N is shown in FigO Since S/L is plot- 
ted, the large- L result is basically the slope s introduced above and seen to increase 
monotonously from 0.6116 to 1 as 8 varies between and 1. Near (5 = one finds 
a special feature due to an eigenvalue Wk ~ 8 2 of ps which leads to a non-analytic 
dependence Ss(L) ~ S 2 In \S\. 

Finally, we have investigated the transverse Ising model on a ring with 
Hamiltonian 

N 

and calculated the sublattice entanglement in both representations. The result for the 
critical case ft, = 1 is very similar to that for the XX model and shown in Fig[2] on 
the right for up to N — 14 sites. The curves can again be fitted by a linear function 
and give s = 0.258(2) and sq = 0.700(1) in the spin picture while Sp is exactly 
proportional to L with a slope s — 0.6001. Thus both entropies are extensive but the 
fermionic value is more than twice the spin value. 

More interesting is the non-critical case. Fig. 4 shows on the left the results for 
the spin representation. One sees that for the larger sizes, Ss/L shows a a maximum 
near the criticial value h = 1. In the disordered region h > 1, the curves approach a 
limit and S is therefore extensive, but in the ordered region h < 1, they depend on L 
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and S is not extensive. Rather one has Sg — > In 2 for h — > 0. This is the same value as 
for a subsystem in the form of one block and has the same origin. The ground state 
in this limit is a superposition of the two states with all spins having a x — +1 and 
all having a x = —1, respectively. This GHZ state leads to a RDM with two non-zero 
eigenvalues w = 1/2 and thus to In 2 for the entropy. 




Figure 4. Sublattice entanglement entropy per site in TI rings as a function of 
the transverse field. Left: calculated in the spin basis for rings of different size 
N. Right: calculated in the fermion basis in the large- TV limit. Note the different 
vertical scales. 



The fcrmionic case is different and can be treated analytically although the 
ground state is more complicated than in the XX model. Following [T2], H] , one has to 
find the eigenvalues (2£ — l) 2 of the matrix 

M=(2C-1-2F)(2C-1 + 2P) (27) 

where Fjj = (c^cj) and all sites are on one sublattice. However, due to the 
translational invariance of the subsystems, one can work in momentum space. Then 
the matrix M decomposes into 2x2 blocks involving the correlation functions 



C q — (a q a q ), C- q — (aL g a_ g ) and F q = (a q a_ q ). As a result, one finds 



1 



1 



[(1 + h 2 ) 2 - 4/j 2 cos 2 ?] 1 ^ 



from which the entanglement entropy is obtained as 

*J=- >: [C„lnC 9 +(l-C 9 )Ml-C,)] 



E 

0<q<7r/2 



(28) 



(29) 



This is always extensive and shown in Fig. 4 on the right. For h > 2.5 the curve 
practically coincides with the one in the spin representation. This can be attributed 
to the short correlation length whereby only neighbouring sites see each other. The 
asymptotic form is S/L ~ In h/h 2 since C, q = 1 — 0(1 /h 2 ) for large h. As h is reduced, 
however, the curve rises continuously to one, corresponding to the value Sf — Lin 2, 
with no sign of the long-range order. For three values of h, the ( q are independent of 
q, namely for h — (( q = 1/2, Sf 
where ( q = (1 + l/y/2)/2 and 



L In 2), for h = oo = 1, S F = 0) and for h = 1 



Sf = L 



In 2- — ln(l + V2) 
V 2 



(30) 
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This gives the value s = 0.6001 cited above. Near h = 1, the entropy per site varies 
asymptotically as (h — 1) In \h— 1| and thus has infinite slope. This is the signature of 
the phase transition in Sp. The interpretation of the extensivity is similar as in the 
XX model, but in the ground state 



0<q<ir 

one has now a coupling of the two single-particle states (q, —q) in the two sublattices 
if one inserts (|24p . The corresponding operators appear twice in the product, since 



Since Tl and XX chains are related by a dual transformation, one could expect 
a relation between the corresponding entanglements. If the subsystem is a block in 
a ring, such a connection indeed exists [57] ■ m the sublattice case, however, we have 
not found a similar result. 

5. Conclusion 

We have studied the entanglement in spin chains for the case that the subsystem is 
not singly connected. We demonstrated that working in the spin and the fermion 
representation leads in general to different RDM's and to different entanglement 
entropies. We did this by looking at two extreme cases, namely a subsystem of only 
two sites and one in the form of a whole sublattice. The first one displays the effect 
particularly clearly, while the second one provides an example, where it is particularly 
large. 

At the level of the wave function, there is no difference between the two repre- 
sentations. One can rewrite the spin expression directly into the occupation-number 
form. However, the operators sample different information. In the fermion picture, the 
spin function (af <rj) needs all sites between i and j, whereas in the spin picture the 

same is true for the fermion function (c\cj). If some of these sites do not belong to the 
chosen subsystem, the reduced density matrices giving these expectation values will 
usually not coincide. If one determines them by integrating out degrees of freedom, 
this difference arises, because one needs Grassmann variables in the fermionic case 
[2"5I [T5] . This can lead to sign changes in the terms contributing to a particular ma- 
trix element, as compared to the spin calculation, and thus to a different final result. 
One can see this explicitly by considering an XX chain with four sites. Thus while the 
eigenvalues of the spin RDM are directly related to the coefficients in the Schmidt de- 
composition of the state, this does not necessarily hold for those of the fermion RDM. 
They and the resulting entropy measure the entanglement in a somewhat different way. 

Note added. The difference between the two representations discussed here has 
also been noted in a preprint by V. Alba et al., arXiv:0910.0706 ( which just appeared. 
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